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ABSTRACT 


Turbulent  pipe  flow  experiments  have  been  obtained  under  conditions 
which  were  similar  to  the  peak  heating  conditions  of  high  performance 
ballistic  missiles  (approximately  20  percent  lower  enthalpy  and  one- 
half  the  stagnation  pressure).  Several  typical  ablation  materials  were 
investigated  to  determine  their  performance  under  these  conditions. 

It  was  possible  to  determine  the  effective  heat  of  ablation  for  each 
of  these  materials  and  to  experimentally  demonstrate  the  difference 
between  the  ablative  process  in  laminar  and  turbulent  flow.  In  this 
paper  the  Teflon  experiments  are  discussed  in  detail  to  demonstrate 
the  validity  and  power  of  this  technique. 
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SYMBOLS 


B  defined  by  Equation  (10) 

Cp  specific  heat 

D  inside  pipe  diameter 

f  fraction  of  ablated  material  ^vhich  vaporizes 

h.  enthalpy  or  energy  per  unit  mass 

h.^  heat  of  vaporization 

(Ah)^  enthalpy  difference  across  boundary  layer 

effective  heat  of  ablation,  excluding  radiation  emission 
«e£f  effective  heat  of  ablation  including  radiation  emission 
L.  length 

PH  ratio  across  boundary  layer 
m  mass 

m  mass  rate 

M  molecular  weight 

n  exponent  defined  by  £q.  (2)  or  £q.  (3) 

q  heat  transfer  rate 

RT^  reference  energy,  33.86  Btu/lb 

R  inside  pipe  radius 

R^  Reynolds  nvunber,  pviD/fi 

t  time 

T  temperature 
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U  gas  velocity 

X  percent  phenolic 

|3  transpiration  factor 

€  emissivity 

H  viscosity 

$  ratio  of  heat  transfer,  with  and  without,  mass  injection 

p  density 

A  Prandtl  number 

SUBSCRIPTS 
av  average 

cw  cold  wail 

f  final 

i  initial 

j  with  injection 

laminar 
m  material 

rad  radiation 

s  stagnation 

t  turbulent 

V  vapo  r 

w  evaluated  at  wall 
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Introduction 


«• 


The  problem  of  dissipating  large  quantities  of  heat  by  ablation  has 
received  considerable  attention  in  the  literature  during  the  past  few 
years.  Ablation  is  presently  the  favored  method  of  handling  the  re-entry 
heating  problem  for  ballistic  missile  nose  cones.  It  is  being  considered 
for  use  in  large  rocket  nozzles  and  is  one  of  the  methods  of  protecting 
satellites  from  aerodynamic  heating  during  re-entry. 

A  considerable  amount  of  theoretical  and  experimental  work  on 
this  subject  has  been  published  ^ The  theoretical  understanding 
of  the  ablation  process  in  pure  materials,  such  as  quartz-like  substances, 
has  reached  a  high  degree  of  sophistication^*^.  Combination  of  theory, 
experimental  and  flight  test  results  have  produced  gratifying  results  7*  o. 
However,  all  this  work  with  the  exception  of  the  flight  tests,  has  been 
restricted  to  the  laminar  ablation  process.  From  a  look  at  the  Rey¬ 
nolds  number-heat  transfer  history  of  a  high  performance  ballistic  mis¬ 
sile  nose  cone,  it  can  be  seen  that  95%  of  the  material  ablated  from 
the  body  is  lost  in  a  turbulent  flow  region. 

The  ablation  process  in  a  turbulent  flow  is  considerably  different 
from  laminar  ablation.  This  is  particularly  true  for  those  materials 
which  depend  on  vaporization  for  obtaining  their  high  values  of  effective 
heat  of  ablation.  (^alitatively,  this  difference  is  the  relative  effective¬ 
ness  of  the  vaporized  material  in  altering  the  boundary  layer  profile  in 
such  a  manner  as  to  reduce  the  temperature  gradient  at  the  surface, 
i.  e.  ,  heat  transfer.  In  a  laminar  flow  the  injected  material  is  very 
effective  in  producing  lower  gradients  at  the  surface.  However,  in 
turbulent  flow,  with  its  fuller  boundary  layer  profile  and  its  turbulent 
motion  closer  to  the  body,  the  injected  material  is  far  less  effective. 

Turbulent  ablation  is  sufficiently  different  so  that  experiments  need 
to  be  performed  in  support  of  the  theoretical  predictions.  To  date  no 
significant  laboratory  experiments  have  been  performed  where  both  the 
Reynolds  number  and  the  gas  enthalpy  were  high  enough  to  produce  tur¬ 
bulent  ablation.  To  fill  this  gap  in  our  knowledge,  the  large  high  pres¬ 
sure  arc  plasma  generator  was  developed.  Still  the  requirements  of 
duplicating  both  the  flight  enthalpy  and  Reynolds  number  of  high  per¬ 
formance  nose  cones  coupled  with  the  requirement  of  geometric  scaling 
required  arc  facilities  beyond  the  present  state  of  the  art.  To  this  end 
large  high  pressure  arc  facilities  have  been  developed 9. 

By  relaxing  the  requirement  of  geometric  scaling,  one  is  able  to 
design  significant  experiments  at  power  levels  which  are  practical  with 
presently  available  plasma  generators  by  utilizing  a  turbulent  pipe  flow 
as  the  test  media.  It  takes  less  power  (and  Reynolds  niunber)  to  pro¬ 
duce  a  turbulent  pipe  flow  than  it  does  to  produce  a  turbulent  boundary 


layer  flow.  The  inside  of  a  pipe,  through  which  the  plasma  is  stream¬ 
ing  at  a  slightly  subsonic  Mach  number^  appeared  like  a  reasonable 
simulation  of  the  flight  situation.  Hypersonic  turbulent  flow  has  been 
shown  to  be  similar  to  incompressible  flow  both  in  the  shock  tube  and 
in  hot  wind  tunnels  ^1.  The  effect  of  pressure  gradient  is  diminished  in 
turbulent  hypersonic  flow  due  to  the  large  density  ratio  across  the  bound¬ 
ary  layer.  The  flat  plate  incompressible  flow  is  a  good  simulation  of  the 
hypersonic  flow  situation.  The  inside  of  the  pipe  or  pipe  flow  is  most 
certainly  a  good  approximation  to  the  flat  plate  case.  Thus  the  pipe  flow 
experiment  was  initiated  with  the  purpose  of  first  establishing  the  validity 
of  and  evaluating  this  technique  for  hypersonic  simulation  and  to  produce 
meaningful  results  on  the  process  of  turbulent  ablation.  This  paper 
describes  this  experimental  technique,  gives  some  results  and  discusses 
their  significance  in  view  of  the  analytical  understanding  of  ablation. 

Turbulent  Pipe  Experiments 

Consider  the  Subsonic  flow  in  a  pipe  just  aft  of  the  plenum  chamber 
of  a  high  pressure  arc.  A  sonic  exit  nozzle  at  the  aft  end  of  the  pipe, 
in  conjunction  with  the  regulated  air  supply  controls  the  pressure  level 
in  the  pipe.  In  this  pipe  either  fully  developed  turbulent  pipe  flow  (at 
Reynolds  numbers  above  about  10,000)  or  a  turbulent  boundary  layer  (at 
Reynolds  numbers  above  500,000,  as  long  as  the  boundary  layer  thick¬ 
ness  is  less  than  half  the  pipe  diameter,  6  /D  =  .  500)  can  be  developed. 
The  flow  in  such  a  pipe  is  a  reasonable  representation  of  turbulent  hyper¬ 
sonic  flow  due  to  the  fact  that  in  hypersonic  turbulent  flows  the  incom¬ 
pressible  flat  plate  relationships  appear  to  be  applicable  10.  Of  course, 
pipe  flow  is  only  a  special  case  of  flat  plate  flows. 

The  difficulty  with  pipe  flow  ablation  experiments  is  the  changing 
pipe  specimen  diameter  with  time  due  to  ablation  and  the  resulting 
changes  in  heat  transfer.  Some  rough  limits  can  be  established  for  this 
type  of  experiment  due  to  this  effect.  The  energy  flux  in  the  pipe  is 
the  mass  flow,  puA  ,  times  the  enthalpy,  hs,  or 

E  =  (/OaA)hs 


The  heat  of  ablation,  is  defined  as 

dR/it 


(1) 


(2) 
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The  turb 


where  is  the  density  of  the  material  to  be  ablated, 
ulent  heat  transfer  rate  was  given  in  Ref.  10  as 

<l-.oa9^R;®A-^"AK 


If  you  are  willing  to  tolerate  heat  transfer  rate  changes  as  great 
as  a  factor  of  two  during  a  given  experiment,  i.  e.  ,  qfinal  ~  Original* 
then  the  diameter  change  of  the  model  specimen  in  terms  of  Df/lJo  must 
be  no  larger  than  1.47.  A  limitation  for  pipe  flow  experiments  can  be 
established  from  a  combination  of  the  above  equations.  To  obtain  quanti¬ 
tative  limits,  values  of  the  density  and  heat  of  ablation  of  the  model 
material  must  be  assumed.  For  a  material  with  a  specific  gravity  of 
two,  a  heat  of  ablation  of  5000  Btu/lb,  and  a  testing  time  of  10  seconds, 
the  heat  transfer  rate  which  will  give  the  limiting  diameter  change  in 
10  seconds  can  be  plotted  against  energy  invested  in  the  air.  Figure  1 
was  drawn  for  a  pipe  with  a  flow  Mach  number  of  0.  5.  It  can  be  seen 
that  somewhat  marginal  experiments  can  be  performed  in  Model  B  (i.  e.  , 
less  than  one  megawatt  of  power  in  the  gas).  In  the  multiple  arc 
generator,  pipe  experiments  can  be  performed  at  the  realistic  heat  trans¬ 
fer  rates  of  1  -  2  KVf/crn. 


The  foregoing  discussion  assumed  utilizing  the  pipe  as  a  fully 
developed  turbulent  pipe  flow.  To  create  such  a  flow  requires  Reynolds 
numbers  based  on  diameter  of  about  10,000  and  a  pipe  length  of  15  to 
20  diameters.  The  long  pipe  needed  to  produce  the  pipe  flow  extracts 
a  significant  amount  of  energy  out  of  the  gas  and  consequently  the  en¬ 
thalpy  level  is  reduced.  With  a  large,  high  pressure  plasma  generator, 
as  described  in  this  paper,  sufficiently  high  Reynolds  numbers  can  be 
achieved  to  have  actual  transition  and  a  turbulent  boundary  layer  before 
the  boundary  layer  fills  the  pipe. 

Calculation  of  the  conditions  in  pipes  are  supplied  with  hot  plasma 
by  the  arc  generator  described  in  Ref.  9  are  shown  in  Fig.  2.  It  can 
be  seen  that  assuming  transition  occurs  at  or  before  a  Reynolds  number 
of  0.  5  X  10°  (based  on  the  length  from  the  plenum  chamber)  that  a  re¬ 
gion  is  available  in  which  a  turbulent  boundary  layer,  which  has  not  yet 
completely  filled  the  pipe,  exists.  This  point  is  indicated  by  6 /D  =  .500 
or  the  line  labeled  "Fully  developed  pipe  flow.  "  The  conditions  under 
which  the  pipe  diameter  will  change  by  a  factor  of  1.47,  i.  e.  ,  qf/qo  =2.0 
are  also  shown.  It  is  not  considered  to  be  practical  to  operate  above 
this  line. 


The  turbulent  pipe  flow  test  is  not  an  ideal  experiment  for  study¬ 
ing  the  details  of  the  turbulent  ablation  process.  However,  it  does  offer 
a  convenient  method  for  evaluating  ablating  materials  under  turbulent 
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boundary  layer  conditions.  Most  important,  it  can  be  performed  at 
power  levels  which  are  high  but  not,  economically  or  otherwise,  in¬ 
feasible. 

Experimental  Procedure  * 

Air  stabilized  arc  plasma  facilities  have  been  described  in  detail 
in  the  literature  both  from  the  point  of  view  of  operational  capability 
and  as  a  tool  for  ablation  measuring  techniques  9. 12, 13, 14  -phe  plasma 
facility  used  for  the  present  experiments  was  described  in  Ref.  9,  and 
is  schematically  shown  in  Fig.  3.  The  facility  consists  of  a  carbon 
cathode,  a  water-cooled  anode,  settling  or  plenum  chamber,  transition 
section,  pipe  test  section  and  a  sonic  exit  nozzle.  The  arc  burns  in¬ 
side  the  anode  (or  the  first  nozzle)  and  is  stabilized  by  tangential  air 
injection.  The  hot  plasma  passes  from  the  arc  into  the  plenum  chamber 
where  unsteady  effects  are  averaged  out,  into  the  transition  section  and 
through  the  pipe  at  a  slightly  subsonic  Mach  number  and  finally  it  ex¬ 
hausts  to  the  atmosphere  through  the  second  or  sonic  throat. 

The  pipe  section  serves  as  the  test  section  for  the  present  ex¬ 
periments.  In  order  to  determine  the  gas  properties  of  the  plasma  in 
the  pipe  for  a  given  test,  an  energy  balance  for  the  entire  system  is 
maintained.  Water  calorimetry  measurements  are  obtained  for  all  sec¬ 
tions  of  the  arc  prior  to  the  pipe  test  section.  Separate  measurements 
are  performed  to  determine  the  heat  transfer  rate  in  the  pipe  by  replac¬ 
ing  the  ablating  specimen  with  a  copper  water-cooled  section.  At  identical 
conditions  this  calorimeter  section  is  replaced  by  the  ablating  material 
under  consideration. 

During  ablation,  the  inside  diameter  of  the  pipe  increases  with 
time  due  to  the  ablation  and  as  a  result  the  heat  transfer  rate  to  the 
pipe  decreases.  To  check  the  relationship  between  the  heat  transfer 
rate,  pipe  diameter  and  the  relative  size  of  the  pipe  and  the  entry  sec¬ 
tion,  calorimetry  experiments  were  performed.  Pipes  were  sized  to 
have  inside  diameters  less  than,  equal  to,  and  greater  than  the  exit 
diameter  of  the  transition  section.  This  diameter  was  held  constant 
for  all  tests.  It  was  possible  from  calorimetry  measurements  to  cali¬ 
brate  the  facility,  to  establish  the  existence  of  turbulent  liow  and  to 
evaluate  the  effect  of  pipe  diameter  change  on  heat  transfer  rate. 

The  arc  facility  was  powered  from  a  1.  5  megawatt  capacity  bat¬ 
tery  bank.  In  the  present  experiment  the  power  in  the  plasma  at  the 
entrance  of  the  pipe  test  section  was  determined  to  be  of  thf  order  of 
0.  7  megawatt.  It  was  possible  to  test  ablation  specimens  under  condi¬ 
tions  which  were  similar  to  peak  heating  conditions  of  high  performance 
ICBM  (approximately  20  percent  lower  enthalpy  and  one -half  the  stagna¬ 
tion  pressure).  Each  material  was  tested  at  a  given  operating  condition 
repeatedly  for  various  lengths  of  time  to  establish  the  variation  of  pipe 
di^tmeter  with  time. 
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The  materials  tested  during  the  present  investigation  were  typical 
re-entry  ablaters.  In  addition  to  the  typical  heat  shield  materials, 
tests  were  made  using  Teflon  which  is  not  suited  to  this  type  of  re¬ 
entry.  These  tests  were  made  to  study  the  effect  of  turbulent  blowing 
and  to  help  evaluate  the  difference  between  the  performance  of  the  other 
materials  in  laminar  and  turbulent  flow.  The  Teflon  experiments  were 
particularly  useful  as  verifications  of  this  technique  of  ablation  testing 
and  will  be  discussed  in  detail  in  this  paper. 

Data  Reduction 


The  expression  for  the  heat  transfer  rate  for  a  fully  developed 
turbulent  pipe  flow  can  be  written  as 


.  oa5/Ou(AK)^ 

Re  * 


(4) 


where  (Ah)Q  is  the  difference  between  the  stagnation  enthalpy  and  the 
gas  enthalpy  at  the  wall  temperature  and  all  the  properties  such  as 
density  p,  velocity  u,  and  Reynolds  number  Re  are  based  on  average 
properties. 

For  identical  conditions  of  mass  flow  and  power  in  the  gas  the 
effect  of  pipe  diameter  on  the  heat  transfer  rate  is  from  Eq.  (4), 


'li/lf  (Df/Dl)'* 


The  results  of  the  calorimetry  experiments  are  shown  in  Fig.  2. 
Measured  values  of  heat  -transfer  rates  to  pipes  of  three  different  in¬ 
ternal  diameters  are  shown  as  a  function  of  the  power  in  the  air  as 
determined  from  an  energy  balance.  The  experiments  are  performed  at 
constant  air  mass  flow.  The  results  for  the  0.  55  and  0.  75  inch  inside 
diameter  pipes  are  adequately  predicted  by  Eq.  (1).  Furthermore,  the 
variation  of  the  heat  transfer  rate  with  pipe  diameter  is  predicted  by 
Eq.  (5).  This  result  appears  to  be  valid  to  the  point  where  the  test 
section  diameter  is  equal  to  the  transition  section  exit  diameter.  For 
cases  where  the  pipe  section  is  larger  than  the  entry  (a  stepped  entry 
section),  Eq.  (4)  does  not  predict  the  data.  The  step  in  the  entry  sec¬ 
tion  produces  vorticity  which  increases  the  heat  transfer  rate  to  the 
pipe  immediately  downstream  of  the  step  to  values  greater  than  those 
predicted  for  turbulent  pipe  flow,  Eqs.  (4)  and  (5).  The  variation  of 
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heat  transfer  rates  for  a  stepped  configuration  can  be  empirically  rep¬ 
resented  within  the  accuracy  of  the  data  as 


The  agreement  of  the  results  of  the  calorimetry  with  Eq.  (1) 
indicates  that  turbulent  pipe  flow  was  achieved  during  the  present  tests. 
The  transition  section  has  an  L/D  of  8  and  transition  to  turbulent  pipe 
flow  is  apparently  achieved  in  this  section.  Furthermore,  the  impor¬ 
tance  of  establishing  the  variation  of  heat  transfer  rate  with  changes  in 
pipe  diameter  is  demonstrated  by  the  effect  of  the  entrance  step  as 
given  by  Eq.  (6).  The  effect  of  change  in  pipe  diameter  can  be  reduced 
by  the  use  of  a  facility  capable  of  delivering  a  larger  amount  of  power 
to  the  air  and  one  which  consequently  can  accommodate  larger  sample 
sizes.  Such  a  facility  has  been  developed  9  and  allows  turbulent  pipe 
tests  to  be  undertaken  with  small  heat  transfer  variations,  small  cor¬ 
rections,  and  consequently  much  more  accuracy. 

To  establish  the  effective  heat  of  ablation  of  the  various  materials 
about  ten  samples  of  each  material  were  run  at  reproduced  conditions 
(plus  or  minus  5  percent  in  enthalpy)  and  the  weight  of  each  specimen 
was  measured  before  and  after  the  test.  The  weight  of  material  ablated 
during  the  test  is 

-  dm  -  Z'nr  RoiR(OmL  (7) 


where  the  negative  sign  indicates  a  loss  of  material  with  an  increase  in 
pipe  radius.  Integration  of  Eq.  (7)  between  the  initial  and  final  condi¬ 
tions  results  in  the  following  expression 

Rf /Pi  r  I  (8) 

N  TTR-^n>mL 


From  the  measurement  of  the  weight  loss  during  the  test  and  use  of 
Eq.  (8),  the  radius  change  of  the  specimen  was  calculated.  This  value 
represents  an  average  pipe  diameter  of  the  specimen  after  the  test. 
Several  samples  were  sectioned  and  measured  in  detail  to  determine 
the  change  in  pipe  radius  as  a  function  of  axial  distance  along  the  pipe. 
The  values  calculated  from  Eq.  (8)  were  compared  and  found  to  ade¬ 
quately  represent  true  average  pipe  radius.  The  effective  heat  of 
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ablation  without  re  radiation*  can  be  written  as 


(0) 

* 


/Om 


(9) 


Experiments  were  conducted  for  different  exposure  times  so  that 
for  each  of  the  materials  tested  a  curve  of  radius  change  with  time 
could  be  drawn.  Conditions  were  controlled  so  as  to  result  in  repro¬ 
ducibility  from  run  to  run  within  approximately  5  percent.  The  cold 
wall  heat  transfer  rate  to  the  pipe  section  was  determined  from  the 
calibration  curves,  Fig.  2.  Since  there  were  slight  variations  of  power 
in  the  Jir  at  the  test  section,  i.  e.  ,  the  enthalpy,  the  cold  wall  heat 
transfer  rate  values  for  a  given  material  test  sequence  were  averaged. 
The  average  value  of  the  stagnation  enthalpy  was  also  used.  In  order 
to  correct  the  cold  wall  heat  transfer  rate,  qcw>  to  the  wall  tempera¬ 
ture  during  ablation,  the  surface  temperature  of  the  ablating  material 
was  needed.  It  was  not  possible  to  measure  this  temperature  during 
the  present  investigation.  However,  the  vaporization  temperature  of 
quartz  as  a  function  of  pressure^^  indicates  that  at  a  pressure  of  10 
atmospheres  a  wall  temperature  of  6400°R  should  be  expected  for  quartz¬ 
like  materials.  During  ablation,  materials  flow  before  evaporating  and 
tend  to  reduce  this  temperature  somewhat.  Theoretical  calculations  by 
Hidalgo®  indicate  temperatures  of  the  order  of  5800°R  in  the  nose  re¬ 
gion  under  conditions  similar  to  the  present  experiments.  As  a  result, 
the  wall  temperature  for  the  present  tests  was  assumed  to  be  6000®R 
for  all  materials  except  Nylon  Phenolic  and  Teflon,  both  of  which  will 
be  discussed  separately  in  a  later  section  of  this  paper.  This  temper¬ 
ature  was  used  to  perform  the  hot  wall  heat  transfer  correction, 

(hs  -  CpTw)/hs. 

A  typical  plot  showing  the  variation  of  radius  with  time  (curves 
are  for  Teflon  tested  at  two  different  values  of  enthalpy)  is  shown  in 
Fig.  3.  The  curvature  of  the  trend  of  the  data  is  due  to  the  change 
in  radius  and,  hence,  the  heat  transfer  rate  with  time.  The  local  slope 
of  the  curves  and  local  radius  were  determined  for  discrete  time  inter¬ 
vals.  These  two  values  complete  the  data  needed  for  Eq.  (9).  The 
local  slope  is  dR/dt  and  the  local  diameter  change  determines  the  ratio 
of  the  initial  to  instantaneous  diameter  needed  to  correct  the  heat 
transfer  rate  for  diameter  changes.  The  exponent  used  was  determined 
from  either  Eq.  (5)  or  (6),  depending  upon  the  initial  radius  of  the 
specimen.  In  this  way,  the  value  of  the  effective  heat  of  ablation  was 
determined  for  each  of  the  materials  and  at  several  enthalpy  levels  for 
some  materials. 


’S’  The  effects  of  radiation  in  this  experiment  are  discussed  in  the  next  section. 
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Discussion  of  Results 


Before  discussing  the  results  of  these  experiments,  let  us  con¬ 
sider  the  difference  between  the  laminar  and  turbulent  ablation  processes. 
It  is  now  generally  accepted  that  the  valv^  of  a  material  as  an  ablative 
heat  protection  system  is  measured  in  terms  of  its  effective  heat  of  * 
ablation,  i.  e.  ,  the  heat  transfer  or  energy  absorbed  per  unit  mass  of 
the  ablating  material. 

For  the  case  of  a  subliming  material,  one  which  goes  directly 
from  the  solid  to  the  gaseous  state  (or  one  which  forms  a  liquid  all 
of  which  vaporizes)  the  effective  heat  of  ablation  can  be  written  (from 
Adams  ^ ), 


H 


% 


Cp  Tyy  ♦  Ky  *  C^^)o 


(10) 


where  =  €  ATyy^ 


The  first  term  in  Eq.  (10)  is  the  energy  absorbed  by  the  heat 
capacity  of  the  material  as  it  is  heated  to  the  temperature  of  the  solid 
gas  interface,  T^,.  The  second  is  the  heat  of  vaporization  of  the  ma¬ 
terial.  The  third  term  represents  the  effect  of  blowing  by  the  vaporized 
material  as  it  enters  the  boundary  layer.  It  is  the  only  boundary  layer 
parameter  in  Eq.  (10)  and  consequently  the  only  one  which  can  be 
affected  by  transition.  The  transpiration  factor,  p  ,  depends  on  the 
relative  molecular  weights  of  the  vaporized  ablating  material  and  the 
working  medium,  i.e.  ,  air.  Because  it  is  a  boundary  layer  effect,  it 
has  a  different  value  for  laminar  and  turbulent  flow.  The  denominator 
1  ~  <lrad/4  represents  the  effect  of  re-radiation  from  the  hot  ablating 
surface  to  the  gas. 

For  a  non -subliming- material,  i.e.,  one  which  has  a  liquid  phase 
before  it  vaporizes,  Eq.  (10)  can  be  rewritten 


H 


«ff  * 


CpTyy  f  Ky  ♦  (  A  K  )p 


(11) 


where  f  now  represents  the  fraction  of  the  ablated  material  which  vapor¬ 
izes  while  1-f  remains  a  liquid  and  flows  off  the  body. 
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The  difference  between  laminar  and  turbulent  ablation  can  be  seen 
from  a  consideraticm  of  the  transpiration  parameter  p  .  By  considering 
the  data  of  Baron,  ^  '  Pappas,  ■  and  Rubesin,  Hidalgo®  has  established 
correlations  for  the  reduction  of  heat  transfer  due  to  mass  injection.  By 
extending  the  similar  solutions  of  Kemp,  Hidalgo  calculated  the  reduction 
in  the  heat  transfer  for  laminar  flow  to  be 

»  t  -  0.56  6  (12) 


where  qj  and  q  are  the  heat  transfer  rates  with  and  without  injection  , 
respectively,  ^  is  the  y# ratio  across  the  boundary  layer  without 
injection,  A  is  Prandtl  number  and  B  is  the  blowing  parameter  defined 
as 
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,  (13) 


The  above  equation  was  derived  from  numerical  solutions  of  the 
boundary  layer  equations  for  dissociated  air.  There  is  also  an  effect 
due  to  the  molecular  weight  of  the  injected  gases  as  shown  by  the  re¬ 
sults  of  Baron.  Hidalgo's  correlation  showed  that  the  blowing  parameter 
B  should  be  multiplied  by  {29/M)^f^  where  M^i  is  the  molecular  weight 
of  the  injected  material.  Thus  the  transpiration  factor  for  laminar  flow 
becomes 


29/m;^**  (14) 

A  correlation  from  the  data  of  Baron  for  the  case  of  ^  =1.0 

and  no  dissociation  gave 


/0^  .  0.67  (29/m);^* 


(15) 


A  rough  check  on  the  value  of  Pt  under  re-entry  conditions  from 
the  ablation  data  of  Georgiev  ^  using  Teflon  gives  a  value  of  Pu  of 
0.  60  44.  Thus  the  laminar  transpiration  factor  is  pretty 

well  est^lished  both  experimentally  and  theoretically. 

For  mass  injection  into  a  turbulent  boimdary  layer,  Hidalgo  has 
found  the  following  correlation  equation  to  be  in  reasonable  agreement 
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with  the  experimental  data  of  Pappas. 


yr  .  » 1 . 0.2C29/M)*  B 

where  OC  =  0. 33  for  (29/m;^  i  1.0 
=  0.8  for  (29/M;^i  1.0 

°r  .  0.2  (29/M)* 


(16) 


(17) 


For  the  ablating  materials  of  interest  here,  the  0.33  value  of  OC 
is  applicable  because  the  vapor  products  have  a  greater  average  mo¬ 
lecular  weight  than  air. 

The  striking  difference  between  the  laminar  and  turbulent  trans¬ 
piration  effect  is  shown  in  Fig.  4,  where  the  calculated  reduction  in 
heat  transfer  rates  with  blowing  is  plotted  for  the  case  of  Teflon.  Based 
on  the  experimental  data  which  leads  to  the  values  of  oc  indicated 
above,  the  difference  between  laminar  and  turbulent  ablation  would  be 
minimized  if  very  light  average  molecular  weight  vapors  were  formed 
during  the  gasification  process. 

From  the  previous  equations  a  general  conclusion  can  be  drawn 
about  the  choice  of  a  material  for  turbulent  ablation.  Subliming  mate¬ 
rials  or  materials  which  achieve  their  high  Heff  by  having  a  large 
fraction  vaporized  (thereby  depending  strongly  on  blowing)  may  not  be 
very  good  turbulent  ablaters.  A  good  turbulent  ablater  should  get  most 
of  its  effectiveness  from  a  high  heat  of  vaporization,  i.  e.  ,  the  term 

(fhv). 


As  was  seen  from  Eq.  (9),  a  knowledge  of  the  surface  tempera¬ 
ture  during  the  ablation  is  necessary  for  the  analysis  of  the  experi¬ 
mental  data.  As  previously  discussed  the  surface  temperature  deter¬ 
mines  the  value  of  CpTw  and  consequently  the  value  of  the  hot  wall  heat 
transfer  correction  or  the  incident  heat  transfer  on  the  specimen.  The 
heat  transfer  rates  measured  in  the  calorimeter  experiments  must 
always  be  adjusted  before  they  can  be  used  due  to  the  very  high  values 
of  surface  temperature  of  the  ablating  material. 

The  effective  heat  of  ablation  is  heat  transfer  sensitive  because  of 
the  effect  of  re-radiation  of  energy  from  the  hot  material  surface  into 
the  plasma.  For  materials  such  as  the  quartz-like  materials  under 
investigation  the  surface  is  at  a  temperature  of  about  6000^R  at  pres¬ 
sures  of  five  atmosphere  and  higher.  This  represents  a  re-radiation 
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rate  of  the  order  of  400  Btu/ft^  sec  (depending  on  the  emissivity).  Thus 
the  value  of  Heff  would  tend  to  go  to  infinity  as  the  heat  transfer  rate 
to  the  hot  wall  tends  to  go  to  this  value.  This  effect  is  very  apparent 
in  the  case  of  the  cylindrical  surfaces  of  some  re-entry  vehicles.  The 
maximum  heat  transfer  rates  here  are  barely  greater  than  the  re- 
radiation  heat  transfer  rate  due  tr  the  surface  temperature  of  some  of 
the  presently  used  materials.  Thus  the  ablation  rates  measured  here 
have  been  very  small  and  the  resultant  Heff  abnormally  large. 

In  the  turbulent  pipe  experiment  the  effect  of  re -radiation  is 
eliminated  and  the  values  of  are  not  heat  transfer  rate  dependent. 

This  is  not  precisely  correct  because  the  effective  heat  of  ablation  does 
depend  on  the  body  geometry  or  particle  history  in  the  boundary  layer 
to  some  degree.  However,  this  dependence  is  very  weak  for  turbulent 
flow.  This  elimination  of  the  re -radiation  is  caused  by  the  fact  that 
the  gas  is  not  very  opaque  in  the  conditions  used  in  the  experiments^! 
and  consequently  the  re -radiated  energy  goes  through  the  transparent 
gas  and  is  re -absorbed  by  another  portion  of  the  ablating  surface.  The 
absorption  at  the  pipe  surface  is  far  greater  than  the  absorption  in  the 
gas  for  each  path  length.  Thus  the  turbulent  pipe  is  an  experiment  in 
which  the  effect  of  radiation  is  eliminated  and  as  such’  has  added  merit 
and  significance. 

Thus  for  the  turbulent  pipe  experiment  the  effective  heat  of  abla¬ 
tion  can  be  written  as 

(18) 

*Cp  Tw  +  ^hy  + 

From  the  above  equation  it  can  be  seen  that  a  measurement  of 
Hgff  or  the  rate  of  material  ablation  under  known  -heat  transfer  condi¬ 
tions  for  a  material  which  is  known  to  sublime,’  i.  e.  ,  f  =  1,  leaves  the 
transpiration  factor  or  blowing  effectiveness,  P  ,  as  the  only  unknown 
quantity.  This  is  effectively  the  same  procedure  as  has  been  followed 
in  laminar  ablation  experiments. 

Teflon  meets  the  requirement  of  ablating  by  direct  sublimation 
into  a  vapor  and  is  consequently  suitable  for  measuring  the  blowing 
efficiency,  P  .  Laminar  experiments^  yielded  the  result  that 


H,fJ*950  +  .44  (4h), 


(19) 


To  establish  the  value  of  the  turbulent  blowing  efficiency  p  ,  the  value 
of  H^£  was  measured  at  two  different  enthalpies.  These  measurements, 
see  fig.  5,  produced  a  line  with  the  equation 
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f  »(950)  ♦  (0.14)  (4^)0 


(20) 


Analysis  of  the  measured  data  indicates  an  uncertainty  of  plus  (.  05)  for 
the  value  of  the  blowing  efficiency,  p  ,  in  the  above  equation. 

The  line  for  the  above  equation  has  been  drawn  through  the  lower 
enthalpy  data  and  the  known  value  of  the  heat  of  vaporization  of  Teflon^^. 
The  disagreement  between  the  data  points  and  the  curve  can  probably  be 
explained  in  terms  of  the  prediction  of  the  heat  transfer  rate  change 
with  diameter  from  the  calorimetry.  As  can  be  seen  from  Fig.  6,  Teflon 
ablates  very  smoothly  but  the  calorimetry  was  performed  for  a  step 
change  in  diameter.  Thus  the  decrease  in  the  exponent  of  the  heat  trans¬ 
fer  rate  variation  with  diameter  change  from  1.8  to  .300  as  indicated  by 
the  calorimetry  data,  Fig.  2,  is  probably  not  as  strong  for  the  Teflon 
experiment.  This  effect  would  tend  to  lower  the  high  enthalpy  data.  The 
'agreement  of  the  present  result  with  the  value  anticipated  from  the  low 
speed,  cold  air  results  of  Rubesin  is  quite  gratifying  in  view  of  the 
shortcomings  of  the  present  experiments. 

The  agreement  of  the  Teflon  data  with  the  turbulent  blowing  esti¬ 
mates  and  the  close  check  between  the  turbulent  heat  transfer  predictions 
and  calorimeter  data  appears  to  verify  our  contention  that  the  turbulent 
•  pipe  experiment  measures  turbulent  ablation.  Figure  6  shows  a  photo¬ 
graph  of  a  representative  sample  of  Teflon  as  it  appeared  after  the  ex¬ 
periment.  The  samples  were  exposed  for  various  lengths  of  time,  and 
consequently  the  relative  amounts  of  the  apparent  ablation  on  this  picture 
are  not  significant.  Two  samples  of  Teflon  are  shown  to  indiccte  the 
uniformity  of  the  diameter  change  over  a  large  range  of  total  ablation 
distance.  This  uniformity  was  ^not  as  good  for  flowing  materials. 

Some  materials  other  than  Teflon  were  tested  under  this  program. 
The  additional  complexity  of  heat  shield  materials  which  do  not  sublime 
entirely  but  some  fraction  of  whose  products  flow  and  leave  the  sample 
without  fully  evaporating  does  not  effect  this  technique  adversely.  Exam¬ 
ples  of  this  type  of  material  were  refrasil  laminates  filled  with  a  phenolic 
resin  binder.  Such  materials  ablate  very  uniformly  and  leave  a  thin  but 
firm  charred  layer.  There  were  no  signs  of  flowing  material,  although 
analysis  of  the  data  shows  that  only  15%  -  25%  of  the  refrasil  vaporized. 
The  fraction  vaporized  was  calculated  using  a  blowing  factor  from  Eq.  (17) 
and  an  estimated  molecular  we^ht  of  the  vapor  of  about  40  -  50.  The 
latter  was  calculated  by  Shick^^for  the  reaction 

SiO,^SiO  +  '/*  0*  (2» 
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and  can  also  be  estimated  from  the  laminar  data  for  quartz  .  The 
result  of  substituting  this  molecular  weight  into  Eq.  (17)  was  a  blowing 
parameter,  P  of  0.18.  This  value  was  used  for  analyzing  the  perform¬ 
ance  of  refrasil  and  silica  laminate -phenolic  resin  combinations.  The 
fraction  vaporized  for  these  materials  was  calculated  by  substitution 
into 


(O) 


f  =  CpTw  +  xh, 


'(phenolic) 


♦(i-x)fhv,. . +[x+f(i-x)],S(Ah), 


(22) 


'^Csilico) 


Opaque  Quartz 

This  material  was  essentially  pure  hot  pressed  quartz  made  opaque 
by  a  material  suspended  in  it.  The  material  has  a  greyish  appearance 
when  cold.  When  the  material  is  heated  enough  to  flow,  the  impurity 
goes  into  solution  and  combines  with  the  quartz  to  form  a  dark  blue 
color.  The  physical  properties  of  opaque  quartz  are  apparently  similar 
to  those  of  clear  quartz. 

A  characteristic  difference  between  the  ablation  of  this  pure  ma¬ 
terial  and  the  composites  of  silica  and  plastic  was  that  the  flowing  of 
the  liquid  layer  was  very  much  in  evidence  for  opaque  quartz.  Appar¬ 
ently  the  fraction  of  this  material  which  did  not  vaporize,  i.  e.  ,  1  -  f, 
flows  along  the  pipe  surface  and  out  through  the  sonic  nozzle.  This  was 
very  evident  from  a  thick  coating  of  quartz  found  on  the  exit  nozzle  and 
in  the  supersonic  expansion  region  aft  of  the  sonic  nozzle.  This  was 
opposed  to  the  behavior  of  the  other  materials  in  which  the  silica  which 
did  not  vaporize  was  injected  into  air  stream  and  carried  out  through 
the  nozzle. 

The  flowing  quartz  decreased  the  effective  sonic  throat  size  and 
raised  the  pressure  of  the  tests.  The  oqaque  quartz  experiments  were 
consequently  performed  in  a  range  of  pressure  between  12  and  15  atmos¬ 
pheres,  against  an  average  pressure  of  10  atmospheres  for  all  the  other 
materials.  The  models,  when  inspected  after  cooling,  had  a  very  rough 
inside  surface,  particularly  at  the  downstream  end.  A  characteristic 
wave  pattern  was  noticeable.  Whether  this  appearance  is  representative 
of  the  condition  of  the  surface  during  ablation  or  is  caused  upon  cooling 
is  not  known. 

The  performance  of  opaque  quartz  was  estimated  to  Adams^,  as 
shown  in  Fig.  7.  This  estimate  is  based  on  the  results  of  a  large 
number  of  calculations  by  Hidalgo®,  using  the  Bethe-Adams  theory®  for 
different  re-entry  bodies.  A  region  of  zero  pressure  gradient  on  the 
fore  cone  of  these  bodies  is  used  in  the  calculations  to  minimize  the 
effect  of  geometry  of  Heff.  The  line  on  Fig,  7  is  drawn  roughly  through 
the  mean  of  the  values  calculated  in  this  manner. 
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Nylon  Phenolic 


The  ablation  of  nylon  phenolic  is  a  more  complex  process  than  the 
materials  considered  so  far,  involving  not  only  melting  and  evaporation 
but  also  pyrolysis  of  the  plastic  forming  many  different  chemical  species 
and  combustion.  These  processes  finally  leave  a  charred  layer  which 
has  been  completely  stripped  of  its  nylon  and  which  is  extremely  hot  and 
brittle.  This  char  residue  is  due  completely  to  the  decomposition  of 
the  phenolic  resin^. 

The  appearance  of  the  nylon  phenolic  specimen  after  these  experi¬ 
ments  is  quite  characteristic  and  completely  different  from  the  other 
materials  tested.  There  is  no  evidence  of  any  flow  of  nylon  although 
this  has  been  observed  experimentally  in  high  speed  motion  pictures  of 
stagnation  point  ablation  in  a  supersonic,  low  pressure  arc  wind  tunnel. 
The  pipe  had  many  deep  longitudinal  gouges  up  to  l/4  inch  deep.  The 
surface  was  quite  irregular  due  to  these  gouges,  as  can  be  seen  on  the 
photograph,  Fig.  6.  The  charred  layer  was  extremely  brittle  and  in 
many  places  broke  off  during  handling. 

Yet  despite  this  poor  physical  appearance  nylon  phenolic  exhibited 
excellent  effective  heat  of  ablation  values  in  these  tests.  The  analysis 
of  the  nylon  phenolic  ablation  data  was  complicated  by  the  absence  of  a 
reliable  estimate  of  the  surface  temperature  of  the  char  layer  during 
ablation.  From  laminar  experiments*  tlie  effective  heat  of  ablation  was 
found  to  be 


« 1000 +0.52  (Ah), 


Steg^ 


gave  an  estimate  due  to  Scala  which  was 


Hef*)'l500  +  0.50(Ah)o 


(23) 


(24) 


Equations  (23)  and  (24)  indicate  that  the  products  of  the  gasification 
of  the  nylon  and  the  phenolic  resin  have  an  average  molecular  weight  of 
about  50.  In  the  above  mentioned  experiments  the  surface  temperatures 
were  measured  to  be  5000°R.  These  experiments  were  performed  at  a 


*  Unpublished  experiments  by  Georgia v  and  Adams  at  Avco- Everett 
Research  Laboratory. 


fraction  of  an  atmosphere  pressure.  Due  to  the  great  dependence  of 
the  numerical  results  on  the  surface  temperature,  the  interpretation  of 
these  experiments  is  somewhat  speculative  where  surface  temperature 
is  not  well  known. 

One  way  to  explain  the  experimental  values  of  Hgff  for  nylon 
phenolic  is  to  assume  that  the  charred  layer  achieves  a  much  greater 
temperature  than  was  measured,  in  the  laminar  experiment.  Some  spec- 
trographic  measurements  have  indicated  surface  temperatures  greater 
than  5800°R  for  nylon  phenolic  during  ablation.  From  the  ablation 
measurements  H§ff  can  be  calculated  for  any  assumed  surface  temper¬ 
ature  for  the  enthalpy  levels  tested. 

More  accurate  determination  of  the  behavior  of  ablating  nylon  * 
phenolic  will  need  a  good  surface  temperature  measurement  during 
ablation.  However,  from  the  present  experiments  it  can  be  concluded 
that  the  surface  temperature  of  the  char  layer  is  probably  very  high  in 
the  vicinity  of  6000°R.  These  data,  together  with  the  appearance  of 
the  char  layer  after  the  experiment,  lend  emphasis  to  the  important 
position  of  the  char  layer  in  the  behavior  of  nylon  phenolic.  In  order 
to  perform  properly  in  turbulent  flow,  nylon  phenolic  apparently  needs 
a  very  high  surface  temperature.  This  can  only  be  achieved  by  the, 
presence  of  a  complete  and  continuous  char  layer.  If  the  char  layer, 
is  broken  or  removed  even  locally,  then  the  ablation  rate  increases  by 
a  large  factor  until  a  full  char  layer  is  reestablished.  This  may  be  a 
speculative  explanation  of  the  gouges  on  the  inside  of  the  model. 

The  extremely  fragile  state  of  the  char  layer  leads  to  another 
reservation  about  the  use  of  this  data  in  projecting  flight  performance. 
The  present  experiments  were  performed  at  approximately  (within  a 
factor  of  two)  the  proper  pressure,  enthalpy  and  heat  transfer  rate  for 
simulating  ballistic  missile  re-entry.  However,  the  conditions  are 
quite  different  with  respect  to  velocity  and  consequently  also  shear.  The 
shear  on  the  fore  cone  of  the  re-entering  nose  cone  is  approximately 
20  times  higher.  Thus  the  apparent  abi.'ity  of  the  char  layer  to  with¬ 
stand  the  experimental  conditions  of  the  pipe  test  does  not  guarantee 
its  equivalent  behavior  during  re-entry.  Unfortunately  although  experi¬ 
ments  were  performed  at  higher  pressures,  higher  shears  were  not 
possible  in  these  experiments  and  consequently  no  additional  light  was 
shed  in  this  regard. 

Conclusion 


The  data  and  techniques  described  in  this  paper  have  led  to  the 
following  results: 

(1)  Ablation  tests  in  a  turbulent  pipe  flow  have  led  to  a 
measurement  of  the  turbulent  transpiration  parameter. 

(2)  A  niunber  of  interesting  heat  shield  materials  have  been 
evaluated. 
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(3)  The  accuracy  of  the  data  produced  by  the  experiments 
described  is  limited  because  of  limitations  on  presently 
available  arc  plasma  generator.  The  large  corrections 
which  were  necessary  in  the  evaluation  of  the  data  would 
be  minimized  if  a  larger  power  input  were  available. 

The  experiments  described  in  this  paper  have  since  been  con¬ 
firmed  and  extended  in  the  Ten-Megawatt  Arc  Facility^  at  AVCO/RAD^^. 
In  this  facility  the  conditions  shown  in  Fig.  2  have  been  achieved.  That 
turbulent  boundary  layer  flow  has  been  achieved  can  be  shown  by  the 
calorimetric  measurements  shown  in  Fig.  10.  Utilizing  this  facility, 
measurements  were  made  on  Teflon  over  a  large  range  of  conditions. 
These  measurements  confirm  the  relationship  inferred  from  the  sparse 
data  of  this  investigation,  as  shown  in  Fig.  11. 
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ABLATION  EXPERIMENT  IMPRACTICAL 
//  /in  THIS  REGION  /  /  /  / 


Fig-  1  Limitations  of  Turbulent  Pipe  ^Experiment. 
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Fig.  2  Conditions  achievable  in  a  pipe  flow  experiment  with  arc 
delivering  plasma  at  10  atmospheres  and  with  3.  0  MW  of 
power  in  the  air. 
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Fig.  3  Schematic  Drawing  of  Experimental  Equipment. 
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Results  of  Calorimetry  Measurements 
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DATA  THEORY 
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Fig.  7  Teflon  Ablation. 
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Fig.  8  Samples  from  Turbulent  Pipe  Test. 
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